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Abstract 

Shadows  enhance  the  realism  of  computer-synthesized  images.  Highly  realistic  shadows 
can  be  generated  by  ray  tracing,  which  is  known  to  be  very  expensive.  The  Z-buffer  shadow 
algorithm  is  very  attractive  due  to  its  simplicity  and  efficiency,  but  is  restricted  to  point  type 
light  sources  and  hence  is  not  applicable  to  penumbra.  This  paper  presents  a  new  model 
which  approximates  a  simplified  cone  tracing  technique.  The  new  algorithm  utilizes  the  sim- 
plicity and  efficiency  of  the  Z-buffcr  algorithm.  The  errors  introduced  by  the  approximation 
are  analyzed  and  a  solution  to  control  the  errors  is  presented.  Experimental  results  show  that 
the  new  method  is  very  efficient  and  is  suitable  for  many  applications. 


1.  Introduction 

Shadow  casting  techniques  have  received  great  attention  from  computer  graphics  research- 
ers because  shadows  improve  comprehension  and  enhance  the  realism  of  computer- 
synthesized  images.  Shadows  provide  valuable  cues  to  positional  relations  among  light 
sources  and  objects  appearing  in  an  image.  When  used  in  scene  shading,  shadows  greatly 
enhance  information  about  the  shape  of  objects  shown  in  an  image. 

Much  research  has  been  devoted  to  the  stud-  of  shadow  algorithms.  [16]  classified  sha- 
dow algorithms  into  five  basic  types:  Z-buffer,  area  subdivision,  shadow  volume,  preprocess- 
ing, and  ray  tracing.  Much  previous  work  has  been  restricted  to  point  light  sources  [8].  In 
this  context,  shadowing  means  a  classification  of  each  surface  point  as  being  in-shadow  or 
not-in-shadow.  But  in  the  real  world,  a  light  source  is  rarely  a  volumeless  point.  Low  pass 
filter  techniques  [18]  have  been  used  to  soften  the  edge  of  shadows  to  resemble  penumbra, 
but  this  is  far  from  satisfactory  in  most  situations.  More  sophisticated  shadow  algorithms  con- 
sidering both  umbra  and  penumbra  have  been  generated  by  ray  tracing. 
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The  Z-buffer  shadow  algorithm  [22]  has  been  used  very  widely  due  to  its  simplicity  of 
structure,  its  ability  to  comprise  the  visibility  section  of  most  kind  of  Tenderers,  and  its  suita- 
bility for  hardware  implementation.  The  Z-buffer  algorithm  computes  the  light-source  and 
viewpoint  images  by  depth  buffer  algorithms.  Each  pixel  in  the  viewpoint  image  is 
transformed  to  the  light-source  image  and  compared  with  the  depth  in  the  light-source  image. 
U  the  visible  point  at  a  pixel  of  the  viewpoint  image  is  farther  from  the  light  than  the 
corresponding  depth  in  the  light-source  image,  it  is  in  shadow.  The  aliasing  problem  with  the 
Z-buffer  algorithm  has  been  reduced  by  various  sampling  and  filtering  techniques  [18].  But 
since  the  light-source  in  the  Z-buffer  algorithm  is  a  point  source,  it  is  not  applicable  to 
penumbra. 

Ray  tracing  algorithms  generate  excellent  realistic  images,  but  in  general  they  are  compu- 
tationally very  expensive.  Simple  ray  tracing  has  aliasing  problems  as  veil.  [1]  presented  a 
cone  tracing  technique  for  anti-aliasing,  fuzzy  shadows,  and  dull  reflections.  To  improve  the 
efficiency  of  ray  tracing  [15][19]  created  a  hierarchy  of  bounding  vo'umes  to  reduce  the 
number  of  bounding  volume  intersection  checks;  [10][11][14]  employed  3D  space  subdivision 
in  culling  objects;  [12]  [14][20]  made  use  of  the  coherence  which  exists  between  similar  rays. 
Though  a  certain  degree  of  improvement  in  efficiency  has  been  attained  by  the  above  efforts, 
ray  tracing's  computational  expense  still  restricts  its  applications. 

In  this  paper  we  introduce  a  new  approximated  shadow  penumbra  method.  The  new 
penumbra  approximation  model  utilizes  a  modified  Z-buffer  technique,  and  is  quite  simple  in 
structure  and  very  efficient  in  implementation.  Quantitative  analysis  is  given  about  errors 
introduced  by  the  approximation,  so  that  accuracy  can  be  controlled  according  to  application. 

2.  Penumbra  Approximation  Model 

The  intensity  value  at  a  surface  point  in  a  given  environment  viewed  from  a  certain  direc- 
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tion  is  determined  by  emitted  light,  reflected  light,  and  refracted  light  from  that  point.  The 
most  accurate  approach  to  this  problem  known  to  the  authors  is  the  radiosity  method  [6], 
which  accounts  for  all  interreflections  of  light  between  diffuse  and  non-diffuse  surfaces  in  a 
given  environment.  The  high  quality  of  the  images  this  generates  is  paid  for  by  the 
extremely  large  computing  cost  of  both  the  preprocessing  and  the  final  rendering  processes. 
Rather  than  accounting  for  complex  interreflections  within  the  environment  as  does  the 
radiosity  method,  ray  tracing  techniques  consider  primarily  secondary  reflected  and  transmit- 
ted rays  from  mirror  and  refracted  directions.  Although  ray  tracing  techniques  are  not  as 
accurate  as  radiosity  methods,  they  have  been  able  to  produce  strikingly  realistic  images  in 
less  computing  time.  This  paper  is  restricted  to  the  shadow  casting  problem.  We  proceed  by 
first  looking  at  the  approach  taken  to  this  problem  by  previous  ray  tracing  techniques. 

As  shown  in  Fig.  1,  object  0  is  placed  between  the  light  source  L  and  the  surface  S,  gen- 
erating a  shadow  on  the  surface  S.  To  determine  the  illumination  at  point  p  on  the  surface  S 
we  can   place   a  projection   plane   in  space  with   the   point  p   as  the   center  of  perspective 


Figure  1.  A  simplified  cone  tracing  scheme 


Figure  2.  The  projections  of  object  and  light  source 
projection.  Both  object  and  light  source  are  then  projected  perspectively  onto  the  projection 
plane.  Let  Si  and  So  be  the  projected  areas  on  the  projection  plane  of  the  light  source  and  of 
the  object,  respectively  (Fig.  2).  There  are  three  possible  cases.  U  Si  D  So  =  NULL,  then  the 
surface  point  p  is  fully  illuminated.  If  Si  C  So,  i.e.  Si  is  completely  covered  by  So,  then  the 
surface  point  p  is  in  full  umbra.  Finally,  ii  Si  n  So  *  NULL  and  Si  n  So  *  Si,  i.e.  5;  is  par- 
tially covered  by  So,  then  the  surface  point  p  is  in  penumbra.  Thus  the  illumination  level  at 

,       SinSo 
the  surface  pomt  p  is  proportional  to  — - — . 

Si 

In  a  complex  scene,  we  will  surely  have  more  than  one  object  in  the  space.  A  general 
solution  to  compute  the  portion  of  the  light  source  visible  from  any  single  surface  point 
would  be  a  Z-buffer  hidden  surface  removal  algorithm,  using  the  surface  point  as  the  view 
origin.  This  is  a  form  of  cone  tracing. 


Obviously  using  this  method  to  compute  the  penumbra  at  each  visible  surface  point  in  the 
scene  would  be  extremely  expensive.  We  propose  a  similar  but  fair  less  expensive  approxima- 
tion to  this  approach.  The  scheme  shown  in  Fig.  3  is  same  as  the  one  shown  in  Fig.l  except 
that  the  shadow  mask  is  projected  onto  the  projection  plane  by  taking  the  center  of  the  light 
source  as  the  projection  center.  Instead  of  doing  projections  for  objects  in  the  scene  from 
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Figure  3.  The  projection  center  for  objects  is 
the  center  of  the  light  source. 

every  point  on  the  surface  5  at  which  the  shadow  value  is  to  be  computed,  we  propose  doing 
a  single  projection  for  objects  from  the  center  of  light  source  as  in  a  typical  Z-buffer  algo- 
rithm, and  then  using  this  shadow  mask  to  do  cone  tracing  from  the  light  source  for  all  points 
on  S.  This  eliminates  a  large  amount  of  computation.  But  how  accurate  is  the  shadow  gen- 
erated by  this  method?  To  answer  this  question,  we  simplify  the  case  into  two  dimensions  as 
in  Fig. 4.    The  result  we  get  from  this  can  be  directly  applied  to  the  three  dimensional  case. 


In  Fig.  4,  the  light  segment  L  of  length  2*R  casts  the  shadow  of  the  object  O  onto  the  sur- 
face S.  Points  on  S  between  mj  and  uj  are  in  full  umbra.  Points  from  Mj  to  p  i  and  from  u^  to 
P2  are  in  penumbra.  All  points  beyond  pi  ox  pi  axe  fully  illuminated.  A  projection  line  P  is 
placed  in  space  parallel  to  L.  We  take  the  center  of  L,  Oi  as  the  projection  center.  The  seg- 
ment Po  on  P  is  the  projection  of  object  0.  Let  us  look  at  the  lower  portion  of  S.  The  correct 
bound  between  umbra  and  penumbra  should  be  at  uj,  the  intersection  point  of  line  l\  and  S. 
Therefore  the  practical  bound  we  can  get  by  using  mask  /*„  is  mj'.  An  error  A5',  the  shift  of 
the  boundary  between  umbra  and  penumbra,  is  introduced  due  to  taking  the  center  of  the 
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Figure  4.  Errors  introduced  by  this  approximation  model 
light  source  as  the  projection  center  to  generate  the  shadow  mask.  An  important  question  is: 
how  large  is  A5  in  a  given  configuration?  It  can  be  shown  (see  Appendix)  that 


A5  = 


/j*D2*rDi+z)2+D3; 


(1) 


Furthermore,  if  we  have  a  image  plane  /  and  viewing  position  V,  the  final  error  Ai  in  the 


image  caused  by  A5  is 


Aj  =  A5«— ^'cosre;  = = — - — = — 'cosm 


(2) 


For  a  given  environment,  this  configuration  is  fixed.  Therefore,  only  the  projection  plane 
is  movable,  i.e.  only  Dj  is  changeable  in  (1)  and  (2).  Suppose  D2<<D\  and  D2<<03. 
Then  for  any  given  positive  c,  if  we  require  Ai^c,  we  only  need 

n? 


D< 


D4     R*(D^^D-i) 


(3) 


(3)  suggests  that  for  a  given  configuration,  if  D2  satisfies  (3),  the  error  introduced  by  the 
above  scheme  in  the  final  image  is  less  than  the  given  limit  e.  In  other  words,  there  is  always 
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a  range  for  D2  such  that  we  can  maintain  the  results  as  accurately  as  required.  By  doing  a 
similar  analysis  for  the  boundary  between  penumbra  and  fully  illuminated  area  we  can  get 
the  same  result  as  the  one  above.  This  analysis  provides  us  with  an  important  guideline  to 
designing  our  algorithm. 

3.  Algorithm 

We  now  describe  various  different  implementations  of  the  penumbra  approximation 
model  introduced  above.  The  rendering  program  consists  of  two  passes  as  in  Williams's  Z- 
buffer  algorithm.  The  first  pass  creates  the  depth  map  structures  for  each  light  source;  and 
the  second  pass  renders  the  scene,  using  such  depth  map  structures  to  compute  shadows. 

3.1  First  Pass:  Creating  the  Depth  Maps 

Our  major  goal  is  to  make  the  shadows  to  be  generated  as  accurate  as  possible.  But  in 
general  there  is  always  some  tradeoff  between  accuracy  and  efficiency.  We  would  like  to 
develop  different  implementations  with  variant  accuracy  and  efficiency  to  meet  the  require- 
ments of  different  applications.  Here  we  start  with  the  simplest  implementation. 

3.1.1  A  Simple  Implementation 

As  stated  in  the  previous  section,  minimizing  D2  is  the  key  to  maintaining  accuracy.  But 
the  depths  of  objects  in  the  scene  could  vary  over  such  a  large  range  that  no  single  D2  would 
satisfy  the  condition  (3).  In  such  a  case,  if  efficiency  is  critical  and  a  certain  degree  of  degra- 
dation in  accuracy  is  tolerable,  a  naive  choice  is  to  place  the  projection  plane  in  the  middle 
range  of  the  objects'  depth.  This  position  can  be  found  efficiently  by  bounding  the  objects 
with  a  box.  First,  we  find  the  bounding  box  for  each  object;  then  we  find  a  large  box  bound- 
ing all  of  these  objects'  bounding  boxes.  This  large  box  gives  us  the  range  of  depth.  Let  Zmin 
denote  the  depth  value  of  the  front  face  of  the  bounding  box  as  viewed  from  the  light  source, 
^°d  ^max   the  depth  value  of  the  back  face.   We  then  choose   the  position   to   place   the 


projection  plane  parallel  with  the  front  and  back  faces  to  be  CZmin  +  ZmaxJ/^. 

3.1.2  Depth  Map  Array 

If  the  error  introduced  by  the  above  simple  approach  is  intolerable  for  some  application, 
the  scheme  can  be  refined  by  replacing  the  simple  depth  map  by  an  array  of  depth  maps  in 
the  range  from  Zmin  to  Z^ax- 

This  array  is  generated  by  sequential  placement.  We  start  by  placing  the  first  projection 
plane  at  Zmm-  Using  condition  (3),  we  can  find  the  subrange  d^  associated  with  the  first  depth 
map  in  which  the  desired  accuracy  can  be  maintained.  The  next  projection  plane  is  placed  at 
Zmin  +  '^i-  We  repeat  this  displacement  until  we  reach  Zmax- 

We  also  need  to  maintain  a  Z-list,  Z[n],  for  the  searching  of  the  depth  map  later.  The  Z- 
list  contains  a  set  of  z  values  stored  in  increasing  order.  These  z  values  represent  the  depth 
positions  at  which  the  depth  maps  are  located. 

3.1.3  Pyramid  of  Depth  Maps 

If  we  draw  a  cone  from  the  center  of  the  projection  to  include  all  concerned  objects  in 
space  (Fig.  5),  we  find  that  for  a  projection  plane  located  at  a  certain  depth  z,  any  part  of  the 
plane  outside  of  the  cone  is  useless  for  later  computations.  In  other  words,  the  window  for  a 
depth  map  on  a  projection  piane  doesn't  need  to  exceed  this  cone.  The  closer  is  the  projection 
plane  to  the  center  of  projection,  the  smaller  is  the  size  of  the  window  associated  with  the 
depth  map.  Obviously  in  the  depth  map  array  method,  many  depth  maps  have  boundaries 
exceeding  the  cone,  hence  a  lot  of  data  storage  is  wasted.  This  consideration  motivates  us  to 
develop  a  pyramidal  structure,  which  fits  the  cone  natiirally. 

Pyramidal  data  structures  have  been  employed  in  image  processing  for  the  purposes  of 
data  compression,  pattern  recognition,  etc.    However,  this  technique  has  also  been  demon- 


Figure  5.  A  cone  includes  objects  in  the  scene 
strated  to  be  useful  in  computer  graphics  [21],  for  example,  to  do  bandlimited  texture  map- 
ping, etc. 
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Figure  6.  Placements  of  depth  maps  in  pyramidal  structure 
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Let  the  depth  map  at  the  first  level  in  the  pyramid  have  a  size  of  n*m.  At  the  second 
level  we  would  then  have  a  depth  map  of  •^*y,  and  so  on  for  successive  levels.  Higher  lev- 
els can  be  computed  from  lower  levels  by  using  various  kind  of  filtering  techniques.  The  sim- 
plest filter,  and  the  one  we  implemented  in  our  experiments,  is  local  averaging.  The  location 
of  each  level  is  determined  as  the  following  way.  If  the  first  level  is  at  distance  of  z  from  the 
projection  center,  the  second  level  would  be  at  distance  ZqI2,  and  the  nth  level  would  be  at 
distance  20/2""'. 

One  thing  we  must  account  for  is  that  when  multiple  depth  maps  are  used,  there  are 
discontinuities  between  consecutive  levels.  This  can  produce  undersampling  artifacts.  This 
problem  is  overcome  by  interpolating  between  levels  [21].    This  is  described  in  the  next  part. 

3.2  Second  Pass:  Rendering  Scenes 

The  typical  Z-buffer  algorithm  is  used  for  rendering  a  scene.  Each  surface  point  visible 
from  the  viewer  is  transformed  into  light  source  coordinates  first.  The  illumination  value  at 
the  point  is  then  determined  differently  depending  on  which  type  of  depth  map  structure  is 
chosen. 

3.2.1  Indexing  the  Depth  Map 

If  only  a  single  depth  map  is  generated  for  each  light  source,  there  is  no  need  to  search 
for  a  depth  map.  However,  when  the  multiple  depth  map  structure  is  used,  we  must  use  a 
search  algorithm  to  find  the  proper  depth  map(s)  effectively. 

To  find  the  proper  depth  map  in  a  set  of  depth  maps,  we  first  project  perspectively  the 
given  surface  point  p  onto  the  depth  map  located  at  Z„(m,  which  in  our  pyramidal  structure 
will  have  the  highest  resolution,  taking  the  light  as  the  projection  center.  Let  Zp  be  the  depth 
value  of  p,  and  Zp'  be  the  depth  value  in  the  depth  map  at  the  projected  location  of  p.  If 
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Zp-^Zp' ,  i.e.  p  is  closer  to  the  light,  it  is  fully  illuminated.  Otherwise,  Zp'  is  used  for  search- 
ing the  appropriate  depth  map(s). 

Given  Zp',  a  binary  search  in  the  pre-computed  Z-list  is  performed  to  find  the  proper 
depth  map.  If  for  some  i  (OS(<n)  Z[i]  =  Zp' ,  then  i  is  the  index  of  depth  map  we  need. 
Otherwise,  we  must  have  Z[i]<Zp'  and  Z[\+l]<Zp'  for  some  i.  Thus,  i  and  i+1  are  the 
indices  of  the  two  depth  maps  with  which  an  interpolation  will  be  made. 

In  the  pyramid  structure,  the  depth  map  of  the  first  level  is  at  distance  Zq  from  the  light 

zo 
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source,  and  the  ith  level  is  at  Z,  =     ' "    =  •  •  •  =— -.    For  a  given  Z^',  let  i=  log  2  ^   .    If 

Zq 
logl^P   -«  =  0,  then  i  is  the  index  of  the  depth  map  we  need;  otherwise,    [j'J  and    l/J  + 1  are 

needed  to  do  interpolation. 

3.2.2  Computing  Illumination  Values 

First  we  show  how  to  determine  the  illumination  value  with  a  single  depth  map.  The 
interpolation  method  is  then  described.  Assume  the  illumination  value,  I,  is  between  0.0  to 
1.0,  with  1.0  representing  full  illumination  and  0.0  no  illumination.  A  surface  point  at  which 
the  illumination  is  to  be  determined  is  first  projected  on  the  depth  map  perspectively  as 
shown  in  Fig.  3.  If  the  projected  point  is  out  of  the  depth  map  boundary,  then  I  =  1.0  since  it 
is  not  obscured  by  any  objects  from  the  light  source.  Otherwise,  the  light  source  is  perspec- 
tively projected  onto  the  depth  map  by  taking  Zp  as  the  projection  center.  Let  Si  be  the  pro- 
jected area  of  the  light  source  on  the  projection  plane,  and  Si^isibu  be  the  portion  of  5/  which 
is  visible  from  Zp  as  was  computed  by  the  Z-buffer  hidden  surface  removal  algorithm.  Thus 

the  ratio  value  — ^"'       is  the  illumination  value  I. 
Si 

The  above  approach  is  general  for  arbitrarily  shaped  light  sources.  However,  computa- 
tion costs  can  be  improved  greatly  for  a  rectangular  shaped  light  source,  a  very  common  light 
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shape  in  the  real  world.  Here  we  show  a  very  efficient  computation  using  the  summed  area 
table  technique  for  rectangular  light. 

The  summed  area  table  technique  has  been  applied  to  filtering  texture  in  [9].  Here  we 
demonstrate  that  it  also  can  be  used  as  a  very  efficient  way  to  compute  Z/vuit;,.  Let  us  tem- 
porarily replace  the  depth  map  by  a  shadow  mask.  A  shadow  mask  is  a  two  dimensional 
array  each  of  whose  elements  m^  is  between  0.0  and  1.0  in  value,  where  1.0  means  totally 
blocked,  and  0.0  means  not  blocked.  A  value  between  0.0  and  1.0  represents  a  partially 
blocked  value.  An  rectangular  light  projected  onto  the  mask  will  be  a  rectangle  as  shown  in 
Fig.  7.    Obviously, 


I1S1SI2 
yis;su2 


(4) 


=    (X2-Xi)(y2-yi)  -   Jxi,>i,x2,y2 


In  the  scanline  method,  ^xi.yi.zi.yi  '*  computed  by  a  double  loop  for  each  projected  rectan- 
gle. However,  if  we  precompute  a  summed  area  table  of  the  same  dimension  as  the  mask, 


Figure  7.  Using  a  summed  area  table  to  compute  shades 
for  rectangular  shaped  light 


13 


each  element  of  which  is  a,;  =  so,o,i,j,  then 

yiS;Su2 

0sisj:2  0Si£j:2  Osisii  OSiSjri 

0S;S>2  OS;Syi  0S;S.v2  O^j^yi 

We  need  to  compute  this  summed  area  table  only  once.  Then  for  any  given  rectangle  with 
low-left  comer  at  {x\,yx)  and  up-right  at  {xi,yi),  Sx\,y\,xl,y^  can  be  computed  by  just  two  sub- 
tractions and  one  addition  [9].  But  {x\,y\)  and  (xi.yt)  might  not  be  at  exact  grid  points.  In 
this  case,  we  must  compute  the  summed  area  values  for  each  of  the  four  corners  of  our  rec- 
tangle by  bilinear  interpolation.  The  interpolated  values  are  then  applied  to  (5). 

We  have  described  the  shadow  computation  method  for  a  single  depth  map.  When  using 
multiple  depth  map  structures,  suppose  the  optimal  distance  to  index  a  depth  map  is  at  Z^', 
but  there  is  no  one  map  in  our  precomputed  depth  map  structure  satisfying  this  exactly. 
Surely  we  would  not  want  to  use  a  single  nearest  one  to  do  our  job,  since  this  rounding  off 
solution  could  potentially  generate  aliasing  problems  in  the  final  picture.  Interpolation  is 
again  used  to  solve  this  problem.  Let  i  and  i-l- 1  be  indices  of  two  depth  maps  found  by  the 
method  stated  in  3.2.1.  The  distances  of  these  two  depth  maps  from  the  projection  center  are 
Z,  and  Z,  +  i  respectively,  and  we  have  Z,<Zp'<Z,  +  i.  Then  we  first  use  the  method  for  single 
depth  maps  stated  above  to  compute  the  illumination  value  for  each  depth  map.  If  we  have 
two  illumination  values  /,  and  /, +  i  respectively,  the  interpolated  illumination  value  is  then 
given  by  the  formula 

Z|  +  1  ~  Zn  Lin     ~  Z, 

'  =  ^'*  z      -z     ^  ^■^■%      -z  (^) 

4.  Experimental  Results 

Figures   8   to   14  illustrate  some  experimental  results  of  our   penumbra  approximation 
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algorithm.  In  this  first  test  scene  there  is  a  tall  block  casting  a  shadow  onto  the  floor.  The 
light  is  placed  at  the  front-left-top  side  of  the  block.  The  position  of  the  light  has  been  chosen 
so  as  to  make  the  depth  values  of  the  different  portions  of  the  block  vary  over  a  large  range 
as  seen  from  the  light,  so  that  the  characteristics  of  our  algorithm  can  be  demonstrated 
clearly. 

Figure  8  shows  a  picture  rendered  by  using  the  traditional  Z-buffer  algorithm,  where  the 
light  is  point-type.  This  picture  gives  us  a  good  feeling  about  the  geometrical  configuration  of 
the  test  scene,  and  also  can  be  used  to  be  compared  with  other  pictures  of  penumbra  sha- 
dows. 

We  now  replace  this  point  light  source  by  a  square  shaped  light  source.  Fig. 9  to  Fig. 11 
are  generated  by  using  the  same  light  source,  but  with  different  depth  map  structures,  which 
clearly  show  different  shadow  effects.  Fig. 9. a  to  Fig.9.c  all  use  a  single  depth  map,  but  for 
each  image  the  projection  plane  for  computing  the  depth  map  has  been  placed  at  a  different 
distance  from  the  light.  For  Fig. 9. a,  the  projection  plane  of  the  depth  map  has  been  placed  at 
a  position  near  the  bottom  of  the  block.  Since  the  shadow  mask  is  close  to  the  floor  and  far 
from  the  light,  the  penumbra  area  is  the  narrowest  among  the  three  pictures.  For  Fig.9.b, 
the  projection  plane  of  the  depth  map  has  been  placed  at  the  middle  of  the  block;  and  for 
Fig.9.c,  the  projection  plane  of  the  depth  map  has  been  placed  near  the  top  of  the  block, 
which  produces  the  fuzziest  shadow  of  the  three. 

Fig. 10  and  Fig. 11  use  multiple  depth  map  structures.  Fig. 10  is  generated  by  using  a 
pyramidal  structure  of  only  two  levels.  The  values  between  the  two  levels  are  linearly  inter- 
polated. In  this  picture,  the  shadows  of  the  bottom  part  and  the  top  part  of  the  block  are 
close  to  the  shadows  of  corresponding  parts  in  Fig. 9. a  and  Fig.9.c  respectively;  but  the  mid- 
dle part  is  not  very  similar  to  the  corresponding  part  of  Fig.9.b.  Fig. 11  is  generated  by  using 
a  depth  map  array  with  three  levels,  which  gives  us  the  most  realistic  shadow  effect  of  all. 
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Fig. 14  shows  a  synthetic  ball  casting  a  shadow  penumbra  on  a  real  world  scene.  The  back- 
ground was  digitized  from  an  actual  cloth  drapery  by  using  the  NYU  ratio  depth  sensor  [4]. 
Fig. 12  shows  the  depth  map  of  the  background  as  a  gray  image.  In  addition  to  this  depth 
image,  we  took  an  ordinary  intensity  image  of  the  same  scene  (Fig. 13).  The  original  intensity 
values  of  surface  points  in  the  background  were  superimposed  into  the  picture,  and  the  depth 
values  were  used  to  provide  the  3D  coordinates  of  the  surface  points  for  the  shadow  compu- 
tation. 

5.  Conclusion 

We  have  developed  a  very  efficient  modified  Z-buffer  algorithm  for  approximating  sha- 
dow penumbra.  The  new  algorithm  generates  realistic  shadow  penumbra  with  a  performance 
dramatically  better  than  that  of  ray  tracing  based  penumbra  techniques.  The  accuracy  of  the 
shadow  penumbra  generated  by  our  algorithm  has  been  analyzed,  and  can  be  controlled  by 
users  according  to  their  application  needs.  Further  improvement  for  common  rectangular 
shaped  light  sources  was  achieved  by  using  summed  area  table  techniques. 

Appendix:  Error  Analysis  of  the  Penumbra  Approximation  Model 

Refer  Fig. 4  for  the  following  discussion. 

The  actual  boundary  between  umbra  and  penumbra  on  the  surface  5,  is  at  uj.  Since  the 
shadow  mask  on  the  projection  line  P  is  produced  by  taking  the  center  of  light  L  as  the  pro- 
jection center,  the  occlusion  boundary  on  P  is  shifted  by  At/,  which  in  turn  results  in  the 
boundary  between  umbra  and  penumbra  on  surface  S  being  shifted  by  AS  from  ui  to  u^'. 
Since 

=  -^^  (1) 

D1+D2+D2         D1+D2  ^  ' 

then 
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Note  that 

Ad 
Substituting  (3)  into  (2),  we  have 


-i^,    i.e.    Ad  =  R*—^  (3) 

D2         Di'  Di  ^^ 


Let  V  be  the  camera  position,  and  /  be  the  image  plane.    The  error  AS  and  the  error  Ai  in 
the  image  caused  by  AS  have  the  relationship 


At    ^    AS*cos(Q) 
From  (5)  we  get 


(5) 


i,  ,  i,._±.,„,„  .  ^_.,^^J^,.„'^ -,W  (6) 

When  D2«Di  and  D2«D3,  (6)  can  be  simplified  into 

A  J  =  ; •cos(Q)  (7) 

For  a  given  c>0,  if  we  want 


Ai  =  ; •cos(e)  s  € 

Dl*Ds 


we  only  need 


^^  ,  ,,^ £1 (8) 

^  D4     R*(D^  +  Di) 


e" 


05.        o? 


©4     R*(Di  +  D2) 
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Figure  8:  The  shadow  in  this  picture  is 
generated  by  using  the  traditional  Z- 
buffer  algorithm  with  a  point-type 
light  source. 


Figure  9. a 


Figure  9.b 


Figure  9.c 


Figure  9:  These  penumbral  shadows  are  generated  by  using  a  single  depth  map,  where 

the  projection  plane  of  the  depth  map  is  placed  (a)  near  the  bottom  of  the 
block,  (b)  near  the  middle  of  the  block,  and  (c)  near  the  top  of  the  block, 
respectively. 
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Figure   10:  This  penumbral  shadow   is   generated   by   using   a 
pyramidal  depth  map  structure  of  two  levels. 
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Figure   11:  This  penumbral  shadow   is   generated   by   using   a 
pyramidal  depth  map  structure  of  three  levels. 
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Figure  12:  The  depth  map  of  a  cloth 
drapery. 


Figure  13:  The  intensity  image  of  the 
same  cloth  drapery. 
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Figure  14:  A  synthetic  ball  casting  a  shadow  onto  the  cloth  drapery. 
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